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Starting from the orthorhombic magnetically ordered phase, we investigate the effects of uniaxial 
tensile and compressive stresses along a, b, and the diagonal a+b directions in BaFe2As2 and 
CaFe2As2 in the framework of ah initio density functional theory (DFT) calculations. While - 
contrary to the application of hydrostatic or c-axis uniaxial pressure- both systems remain in the 
orthorhombic phase with a pressure-dependent nonzero magnetic moment, we observe an inversion 
of magnetism at a critical pressure happening when the a and b axes approach the tetragonal 
condition. Moreover, this is followed by a jump in orthorhombicity as well as a selective change of 
orbital occupancy. We estimate elastic constants and discuss our results in view of recent reports of 
modified magnetic and structural transitions in BaFe2As2 under externally applied uniaxial stress 
and some proposed phenomeno logical models. 

PACS numbers: 74.70.Xa,61.50.Ks,71.15.Mb,64.70.K- 



INTRODUCTION 

Since the discovery of high-Tc superconductivity in Fe- 
based materials in 2008, pQ an enormous amount of effort 
has been invested to understand the microscopic behavior 
of these materials. Iron pnictides and chalcogenides be- 
come superconductors either by hole- or electron-doping 
the systems, by application of external pressure or by a 
combination of both. In particular, uniaxial pressure is 
currently being intensively discussed as a possible route 
towards modifying the structural, magnetic and even su- 
perconducting properties of these systems. Uniaxial ten- 
sile stress was initially applied to bilayer Fe pnictides 
like BaFe2As2 in order to detwin the crystals. [2-5] Re- 
cent neutron scattering experiments on BaFe2As2 un- 
der compressive stress along the in-plane b direction 
reported a progressive shift to higher temperatures of 
both the tetragonal-to-orthorhombic phase transition as 
well as the magnetic transition. [6] Also, Blomberg et 
al. observed a significant uniaxial structural distortion 
in BaFe2As2 under tensile stress. [7] More recently, it 
was found that epitaxially strained thin films of FeSe on 
a SrTiOs substrate show an increase in critical supercon- 
ducting temperatures up to 65 K, the highest reported T c 
so far. [8 Clearly, crystal lattice strain plays a key role for 
the magnetic, structural and superconducting properties 
in Fe-based superconductors and a better understanding 
of the microscopic origin of such behavior is desirable. 

In this work we analyze by means of density func- 
tional theory (DFT) calculations the effects of uniax- 
ial compressive stress as well as uniaxial tensile stress 
on the magnetic, electronic and structural properties of 
BaFe2As2 and CaFe2As2 at low temperatures. We mea- 
sure stress in terms of equivalent hydrostatic pressure, 
P = Tr(<x)/3, where a is the stress tensor matrix and 
positive and negative pressures correspond to applying 



compressive and tensile stresses respectively. While no 
sign of a true structural or magnetic phase transition is 
observed in the range of pressures between —2 GPa and 
2 GPa, there is an inversion in the magnetization, i.e. 
exchange of ferromagnetic (FM) and antiferromagnetic 
(AFM) directions, at the stresses where a and b reach 
the tetragonal condition a=b. Interestingly, CaFe2As2 
shows these effects at larger pressures than BaFe2As2- 
This has also important consequences on the orbital d xz 
and d yz occupancies and should manifest in a shift of 
the magnetic ordering temperature as we shall discuss 
below. Our results support recent suggestions [9] that 
a shift in the magnetic ordering temperature Tn should 
be observed under application of stress and we provide 
in this context ab initio-derived estimates for the elas- 
tic constants. Also, a sharp orthorhombic to tetragonal 
phase transition is observed for CaFe2As2 when compres- 
sive stress of 7.7 GPa is applied along a+b. 



COMPUTATIONAL METHODS 

Electronic structure calculations were performed 
within DFT with the Vienna ab initio simulations pack- 
age (VASP) [10] with the projector augmented wave 
(PAW) basis [TT] in the generalized gradient approxi- 
mation (GGA). All our structural relaxations were per- 
formed under constant stress using the Fast Inertial Re- 
laxation Engine (FIRE) [I2j [13]. For this purpose, we 
had to modify the algorithm accordingly (see Appendix) . 
Every 10 steps, we cycled through non-magnetic, fer- 
romagnetic, antiferromagnetic-checkerboard, stripe-type 
antiferromagnetic (along unit cell axis a) and stripe-type 
antiferromagnetic (along unit cell axis b) spin configu- 
rations, and then we continued the relaxation with the 
lowest energy spin configuration. As a final converged 
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magnetic configuration in the orthorhombic phase we had 
always found stripe-type antiferromagnetic order either 
along a or along b, as discussed below. 

The energy cutoff in the calculations was set to 300 eV 
and a Monkhorst-Pack uniform grid of (6 x 6 x 6) points 
was used for the integration of the irreducible Brillouin 
zone (BZ). 

RESULTS AND DISCUSSION 

BaFe 2 As 2 

Starting from the low-temperature orthorhombic 
structure with stripe magnetic order, we performed struc- 
ture relaxations under applied uniaxial tensile and com- 
pressive stresses along a (AFM direction), b (FM direc- 
tion) and the plane-diagonal a+b for both BaFe2As2 and 
CaFe2As2- We simulated pressures in the range between 
—3 GPa and 3 GPa. In the tensile stress range, below 
-2.7 GPa we observe in both systems that for stress along 
a, a sudden expansion in a and contraction in b and c axes 
occurs, and analogously (with a change of axes) a similar 
situation arises when pulling apart along b. This feature 
signals the extreme case of absence of bonding within the 
material, and this pressure range will be excluded from 
further discussion. 

Fig. [T] shows the evolution of lattice parameters for 
BaFe2As2 as a function of uniaxial stress along a, b and 
a+b. At P = GPa a (AFM direction) > b (FM di- 
rection). BaFe2As2 remains in the orthorhombic phase 
with nonzero increasing magnetic moment for large ten- 
sile stress (negative pressure). Pulling apart along the 
(longer) AFM direction a (Fig. [I] (a)) the system ex- 
pands along a, strongly compresses along c and shows 
almost no changes along b; similarly, pulling apart along 
the (shorter) FM direction b (Fig. [T] (b)) b expands, c 
compresses and a shows almost no changes except for a 
pressure P = —0.22 GPa (Fig. [I] (e)). In this case, a « 
b and BaFe2As2 shows a sudden jump in orthorhombic- 
ity where a becomes the shorter axis and b becomes the 
longer axis. This interchange happens with a rotation of 
the magnetic order by 90 degrees. Now, the FM direction 
is along the a axis and the AFM direction is along the b 
axis. We will discuss this feature below. 

Tensile stress along a+b acts similarly on both a and b 
directions, which expand, while the c direction strongly 
compresses (Fig. [l](c) and (f)). 

Under application of compressive stress (positive pres- 
sure), we observe in all three cases a strong expansion 
along c, a compression along the direction of applied 
stress; a, b or a+b, and almost no changes or a slight ex- 
pansion along b and a for the first two cases, respectively. 
A slight expansion along a has also been measured upon 
application of uniaxial stress along b by Dhital et al. [6] . 
BaFe2 As2 remains in all cases orthorhombic with nonzero 



decreasing magnetic moment. Since a > b at zero stress, 
we observe the inversion of axes followed by a jump in 
orthorhombicity and a 90 degree rotation of the magne- 
tization when stress is applied along a at P — 0.22 GPa 
(Fig. [I] (a) and (d)). This inversion of axes, with b > a 
for all higher compressive stress values means that the 
spin configuration shown in the inset of Fig. [I] (a) should 
now be turned by 90 degrees, with b pointing along the 
AFM direction. Such an inversion is also observed for 
compressive stress along a+b at much larger pressures 
of P = 2 GPa. 

Fig. [2] shows the evolution of magnetic moment, vol- 
ume and As height in BaFe2As2 as a function of stress. 
The three quantities show a clearly monotonic behav- 
ior independent of the applied stress direction except for 
small jumps at the pressures P — —0.22 GPa (for stress 
along a) and P = 0.22 GPa (for stress along b) where 
the tetragonal condition is almost fulfilled (a « b) (Fig. [2] 
(b), (d), (f)). Even though DFT calculations overesti- 
mate the value of the ordered Fe magnetic moment at 
P = GPa, it is to be expected that the relative changes 
in magnetic moment should provide a reliable descrip- 
tion of the situation of BaFe2As2 under pressure effects 
as shown in previous studies [T3HT6] . 

Except for the pressures P = —0.22 GPa (for stress 
along a) and P = 0.22 GPa (for stress along b), stress 
always enforces a certain degree of orthorhombicity and 
the system remains magnetically ordered with a decreas- 
ing ordered moment as a function of compressive stress 
(Fig. [2] (a)). Moreover, since the c axis continually ex- 
pands from negative to positive pressures, Ha s increases 
accordingly as a function of stress (Fig.[2](e)). These fea- 
tures have a direct consequence on the electronic proper- 
ties of the system. 

As an illustration, we show the evolution of the Fermi 
surface of BaFe2As2 under application of uniaxial stress 
applied along a in Fig. [3] in the lFe/unit cell equivalent 
Brillouin zone. At ambient pressure conditions (Fig. [3] 
(b)) we observe three hole pockets centered around T and 
M and two electron pockets centered around X. Appli- 
cation of modest tensile stress of 0.07 GPa leads to the 
disappearance of the 3d xy hole pocket around f in the 
k z = plane (Fig. [3] (a) and (b)). On the other hand, 
when compressive stress is applied, the 3d xy hole pocket 
around T significantly increases in size, and additionally, 
at a pressure of 1.7 GPa, small electron pockets, of ma- 
jority 3d xy and 3d z 2 character, appear along the f — M 
directions of the BZ (Fig.[3](c) and (e)). The increase of 
the 3d xy hole pocket size with increasing uniaxial stress 
can be explained by the reduction of Fe-Fe distance along 
the a axis, leading to an increased contribution of 3d xy 
- 3d xy bonding. In fact, the effects of tensile and com- 
pressive stress on the electronic structure shown for the 
example of stress along a can be seen also in our calcu- 
lations for both stress along b and along a+b. 

A more detailed information of the orbital occupation 
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FIG. 1: (Color online) Evolution of the unit cell parameters in BaFe2As2 under application of uniaxial stress in the equivalent 
hydrostatic pressure range [—2 GPa, 2 GPa] (a) along a, (b) along b and (c) along a+b. Panels (d)-(f) show a corresponding 
blow-up of the pressure dependence of the lattice parameters in the range [—0.3 GPa, 0.3 GPa]. Negative pressures correspond 
to tensile stress while positive pressures correspond to compressive stress. Note, that the relationship between axes and iron 
moments shown in the inset in (a) is valid for P || a < 0.22 GPa, in (b) for P || b > —0.22 GPa. For a discussion of the reversal 
of AFM order, see the text. 



can be obtained by analyzing the density of states at the 
Fermi level N(E-p). Applying stress both along a and b 
has the same effect on the total density of states of both 
BaFe2As2 and CaFe2As2, but there is a selective orbital 
order as shown in Fig. [4j N(E-p) is predominantly of 3d xz 
character when a > b and of 3d yz character when a < b. 
This means that the dominant character switches from 
3d yz to 3d xz at & || a « 0.22 GPa, and from 3d xz to 3d yz 
at & || b « —0.22 GPa as expected. 

Clearly, approaching or moving away from the tetrag- 
onal situation a=b by applying external stress will have 
important consequences on the magnetic phase transi- 
tion temperature Tn. Tn will shift to higher tempera- 
tures the further away the system is from the tetragonal 
situation. If we consider, f.i., application of tensile stress 
along a (longer direction), the degree of orthorhombic- 
ity increases and Tn should increase. On the contrary, 
by applying compressive stress along a up to pressures 
where a ~ 6, we expect a shift of Tn to lower tempera- 
tures. This is in agreement with a recent proposal that 
predicts a shift in Tn as a function of stress based on a 
phenomenological model [9j H7] . 



CaFe 2 As 2 

The lattice parameters of CaFe2As2 under application 
of compressive stress along a, b and a+b directions show 



a similar overall behavior compared to BaFe2As2 (see 
Fig. |5| except for some important shifts of the pressures 
at which the system exchanges the FM by the AFM di- 
rections. When stress is applied along the a direction, 
we observe at P = 0.67 GPa a compression of the unit 
cell to the tetragonal symmetry. However, analogously 
to the case of BaFe2As2, this is not followed by the 
suppression of the magnetic moments of iron. In fact, 
the c axis expands through the stress application and at 
P = 0.67 GPa the c lattice parameter in CaFe2As2 is too 
large for the formation of an interlayer As-As covalent 
bond, necessary for a transition to a collapsed tetragonal 
phase and suppression of magnetic moments as observed 
under hydrostatic or c-axis uniaxial pressure. [T3UT5| [T8] 

For tensile stress along the (shorter) b direction, the 
inversion of axes happens at P = —0.33 GPa followed by 
a jump in the orthorhombicity. 

In order to investigate the possibility of a structural 
and/or magnetic phase transition at higher pressures, we 
concentrate now on compressive stress along the diag- 
onal of the a6-plane. We find that orthorhombicity is 
preserved up to 7.7 GPa, where a sharp transition to a 
tetragonal phase appears. This transition is of first-order 
type like the orthorhombic to collapsed tetragonal phase 
transition under application of hydrostatic or uniaxial 
pressure along the c axis [13] but in this case, changes of 
magnetic and structural properties take opposite direc- 
tions; the c axis undergoes a sudden expansion of about 
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FIG. 2: (Color online) (a) Evolution of the magnetic moments 
of iron, (c) of the unit cell volume and (e) of the pnictogen 
height under uniaxial pressure in the range [— 2 GPa, 2 GPa]. 
Panels (b), (d) and (e) show the corresponding blow-up of 
the pressure dependence of these quantities in the range 
[—0.3 GPa, 0.3 GPa]. Negative pressures correspond to ten- 
sile stress while positive pressures correspond to compressive 
stress. 



Sij = [C -1 ]ij by performing linear fits to Ui((jj) and C 
is then obtained by inverting the resulting matrix. For 
BaFe2As2, the elastic constant matrix is 



C 



95.2 ±4.3 

27.3 ±4.8 
43.7 ±4.5 



20.4 ±3.4 
130.8 ±6.1 
47.7 ±4.6 



40.8 ±4.5 
64.0 ± 7.0 
81.0 ±5.6 



GPa (1) 



Utilizing Voigt and Reuss averages [19], defined as 
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it is possible to estimate the bulk modulus. Voigt and 
Reuss averages yield 61.9 ±5.1 GPa and 69.3 ± 7.5 GPa, 
respectively, which is in good agreement with our pre- 
vious estimate [13] and the experimental value of 59 ± 
2 GPa [2D]. For CaFe2As2, the elastic constant matrix is 
given by 



C 



148.7 ±18.5 
63.9 ±21.4 
61.4 ± 14.7 



45.6 ± 12.3 
182.4 ± 18.4 
63.1 ± 11.4 



55.5 ±12.7 
81.2 ±17.5 
68.8 ±11.3 



GPa 



(2) 

which results in bulk modulus of 84.3 ± 14.8 GPa and 
77.7 ± 17.2 GPa using Voigt and Reuss averages, respec- 
tively. Both values are in good agreement with experi- 
mentally determined values of 82.9 ± 1.4 GPa [2T] and 
the estimate based on fits to Birch- Murnaghan equation 
of state [15] . 



9.5%, and a and b axes contract while the iron magnetic 
moments order ferromagneticaly and sharply increase in 
value by around 25%. Interestingly though, contrary to 
the application of hydrostatic and uniaxial pressure along 
c axis, the volume change here is significantly smaller, 
namely an expansion by about 0.9%. 

ELASTIC PROPERTIES 

Utilizing data for the response to the uniaxial stress 
along a, b and c [13] axes we can directly evaluate elas- 
tic constants Cij in BaFe2As2 and CaFe2As2 correspond- 
ing to the orthorhombic deformations. We define elastic 
constants to be such that 

3 

where oi and Uj are stress and strain tensor components 
respectively, and indices i and j can be xx, yy, zz. Strains 
are defined to be u xx = (a — ao)/ao, u yy — (b — bo)/bo and 
u zz — (c— co)/co, where ao, bo and Co are equilibrium unit 
cell dimensions (P = GPa). We first directly obtain 



CONCLUSIONS 

In summary, we analyzed the effects of tensile and 
compressive stress along a, b and a±b on BaFe2As2 
and CaFe2 As2 by means of DFT calculations under con- 
stant stress conditions with the help of the FIRE al- 
gorithm. Starting from the low-temperature magneti- 
cally ordered orthorhombic phase, we find in the pres- 
sure range between —2 GPa and 2 GPa no real struc- 
tural phase transitions in both systems except for a pro- 
nounced orthorhombicity jump in the region where a « 
b accompanied by a 90 degree rotation of the magnetic 
order. FM and AFM directions are interchanged, as are 
the orbital occupations d xz and d yz . This inversion of 
axes is purely of magnetoelastic origin. This phenomenon 
should be observable experimentally. Moreover, we also 
find a slight expansion along the a axis when compres- 
sive stress is applied along b in BaFe2As2, as reported 
recently [6]. Our results agree with recent phenomeno- 
logical theories predicting a shift in Tn depending on the 
stress conditions. Finally we provide estimates for the 
elastic constants of BaFe2As2 and CaFe2As2- The Bulk 
moduli calculated with these estimates are in good agree- 
ment with the experimental measurements. 
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FIG. 3: (Color online) Evolution of the Fermi surface of BaFe2As2 under uniaxial stress applied along the a axis shown in the 
lFe/unit cell equivalent BZ (see [22] for the BZ path definition). Panels (a), (b) and (c) show k z = cuts of the Fermi surface 
at pressures of -0.07 GPa, GPa and 1.7 GPa respectively, while panels (d) and (e) show vertical cuts along the diagonal of 
the BZ for pressures of -0.07 GPa and 1.7 GPa. Grey lines on panels (a), (b) and (c) denote boundaries of the 2 Fe/unit cell 
BZ. 



APPENDIX A 

Within the FIRE [12] algorithm, the energy minimiza- 
tion is achieved by moving the system's position r ac- 
cording to the following equation of motion: 



v(t) 



F(t) 



m 



- 7 (t)\v(t)\[e v{t) -e F{t) ] , (3) 



where e x denotes the unit vector along a?, with x = 
v(t),F(t), t is time and j(t) is a time-dependent fric- 
tion parameter which assures that the system is moving 
down the energy hypersurface in an optimal manner as 
long as the power P(t) = F(t)-v(t) is positive. If P(t) be- 
comes negative, the procedure is stopped and relaxation 
is reinitialized in the direction of the steepest descent. 

In order to use FIRE for relaxation of periodic systems, 
we change from the configuration space of 37V atomic 
Cartesian coordinates r a , a = 1...N, to an extended 
system of 3N + 9 coordinates r a = (s a ,h), consisting 
of lattice vectors which are contained in columns of the 
3x3 matrix h = (hij), and fractional atomic positions 
s a within the unit cell, where Cartesian and fractional 
positions are related by r a = hs a . 

When (|3) is rewritten in terms of coordinates one 
just needs to find the appropriate expression for forces 
F a , that is, derivatives of energy with respect to the 
coordinates r a . Since stress and strain tensors, a and u, 



can be defined through 



a = 



1 dE 



H 



= (l + u)h, 



where I is the identity matrix and H is the lattice matrix 
after an infinitesimal deformation, it is easy to show that 



dE 



E 



dE dr? 



dE 



dr? ds<* 

p 

dE _ dE du 
dh ~ du dh 

so that forces in the extended coordinates are given by 



dE 
d 



a( h ) ="^( h ) ■ 



(4) 



pa 



^ a h ,-Va(h 



(5) 



Forces F a and stresses a are obtained from the electronic 
structure code, and are inserted directly into Eq. ([5|, 
taking into account that a = a ext — & %nt , that is, total 
stress is the sum of internal and external stresses applied 
to the surface of the unit cell. 
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FIG. 4: (Color online) Evolution of the orbital resolved den- 
sity of states of BaFe2As2 at the Fermi level N(Ef) with 
stress (a) applied along a, and (b) along b. Lines joining the 
calculated points are a guide for the eye. 
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FIG. 5: (Color online) Evolution of the unit cell parameters in CaFe2As2 under the application of uniaxial stress in the range 
[—2 GPa, 2 GPa] (a) along a, (b) along b and (c) along a+b. Negative pressures correspond to tensile stress while positive 
pressures correspond to compressive stress. Note, that the relationship between axes and iron moments shown in the inset in 
(a) is valid for P || a < 1 GPa, in (b) for P || b > -0.6 GPa. 



